Crossing the phantom divide with a classical Dirac field 
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Abstract In this paper wc consider a spatially flat 
Friedmann- Robertson- Walker (FRW) cosmological modcl| 
with cosmological constant, containing a stiff fluid and 
a classical Dirac field. The proposed cosmological sce- 
nario describes the evolution of effective dark mat- 
ter and dark energy components reproducing, with 
the help of that effective multifluid configuration, the 
quintessential behavior. Wc find the value of the 
scale factor where the effective dark energy component 
crosses the phantom divide. The model we introduce, 
which can be considered as a modified ACDM one, is 
characterized by a set of parameters which may be con- 
strained by the astrophysical observations available up 
to date. 

Keywords Classical Dirac fields, phantom divide 



1 Introduction 

According to the standard cosmology the total energy 
density of the Universe is dominated today by both 
dark matter and dark energy densities. The dark mat- 
ter, which includes all components with nearly vanish- 
ing pressure, has an attractive gravitational effect like 
usual pressureless matter and neither absorbs nor emits 
radiation. The dark energy component in general is 
considered kind of vacuum energy with negative 
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pressure and is homogeneously distributed and, unlike 
dark matter, is not concentrated in the galactic halos 
nor in the clusters of galaxies. The observational data 
provide compelling evidence for the existence of dark 
energy which dominates the present-day Universe and 
accelerates its expansion. 

In principle, any matter content which violates the 
strong energy condition and possesses a positive energy 
density and negative pressure, may cause the dark en- 
ergy effect of repulsive gravitation. So the main prob- 
lem of modern cosmology is to identify the form of dark 
energy that dominates the Universe today. 

In the literature the most popular candidates are cos- 
mological constant A, quintessence and phantom mat- 
ter. Their equation of state is given by w — p/p, where 
w = — 1, w > —1 and w < —1 respectively. Dark en- 
ergy composed of just a cosmological term A is fully 
consistent with existing observational data. However, 
these data do not exclude the possibility of explain- 
ing the observed acceleration with the help of phantom 
matter. The cosmological constant can be associated 
with a time independent dark energy density; the en- 
ergy density of quintessence scales down with the cos- 
mic expansion, and the energy density of phantom mat- 
ter increases with the expansion of the Universe. 

Mostly, the attention has been paid to dark en- 
ergy as high energy scalar fields, characterized by a 
time varying equation of state, for which the potential 
of the scalar field plays an important role. Among 
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collisionless particles (IWanget al. 1 120051 : iMuller 1 12005 ; 
iRabev I [20051 : ICopeland et al. II2006I ). we have the pos- 
sibility of giving to cold dark matter content an origin 
based on the nature of the CDF. On the other hand, the 
stiff fluid is an important component because, at early 
times, it could describe the shear dominated phase of 
a possible initial anisotropic scenario, dominating the 
remaining components of the model. 

The organization of the paper is as follows: In Sec. 
II we present the dynamical field equations for a FRW 
cosmological model with a matter source composed of 
a stiff fluid and a CDF. In Sec. Ill the behavior of 
the dark energy component is studied. In Sec. IV we 
conclude with some remarks. 
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2 Dynam ical field equations 

plored the so called kinetic k-essence models (IChimento and Feinste in | 
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20041 ) . i.e. cosmological models with several k- fields in 



We shall adopt a spatially flat, homogeneous and 
isotropic spacctimc described by the FRW metric 



which the Lagrangian does not depend on the fields 
themselves but only on their derivatives. It was shown 
that the dark energy equation of state transits from a 
conventional to a phantom type matter. Note that for- 
mally, one can get the phantom matter with the help of 
a scalar field by switching the sign of kinetic energy o f 
the standard scalar field Lagrangian (ICaldwellll2002h . 
So that the energy density p p h — — (l/2)$ 2 + V(&) and 
the pressure p p h ~ -(l/2)$ 2 - V($) of the phantom 
field leads to p v ] x + p p h = — $ 2 < 0, violating the weak 
energy condition. 

In the Universe nearly 70% of the energy is in the 
form of dark energy. Baryonic matter amounts to only 
3 — 4%, while the rest of the matter (27 %) is be- 
lieved to be in the form of a non-luminous component of 
non-baryonic nature with a dust-like equation of state 
(w = 0) known as cold dark matter (CDM). In this case, 
if the dark energy is composed just by a cosmological 
constant, then this scenario is called ACDM model. 

Below, we analyze a FRW universe having cosmo- 
logical constant and filled with a stiff fluid and a clas- 
sical Dirac field (CDF). With this matter configura- 
tion, we will sec that the FRW universe evolves from a 
non-accelerated stage at early times to an accelerated 
scenario at late times recovering the standard ACDM 
cosmology. The CDF may be justified by an impor- 
tant property: in a spatially flat homogeneous and 
isotropic FRW spacetime it behaves as a "perfect fluid" 
with an energy density, not necessarily positive definite. 
This pressureless "perfect fluid" can be seen as a kind 



ds z = dt z 



l (t) (dx 2 + dy 2 + dz 2 ) , 



(1) 



where a(t) is the scale factor. The spacetime contains 
a cosmic fluid composed by (i) a stiff fluid p s = p s o/a e 
and (ii) a homogeneous classical Dirac field ip. The 
Einstein-Dirac equations are 



SH 2 -A=^f + p D , (2) 
a 

(r*V< - a) V - 0, (3) 

where H = a/ a is the Hubble expansion rate, a is a 
constant and the dot denote differentiation with respect 
to the cosmological time. Here po represents the energy 
density of the CDF. 

The dynamical equation for the CDF in curved 
spacetime can be obtained using the vierbein formal- 
ism. So, T l are the generalized Dirac matrices, which 
satisfy the anticommutation relations 



{r\T fc } = -2g lk I 



(4) 



with the metric tensor g lk and / the identity 4x4 ma- 
trix. They can be defined in terms of the usual repre- 
sentation of the flat space-time constant Dirac matrices 

7* as 
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where the Dirac matrices 7* can be written with the 
Pauli matrices cr' 3 as 
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with / the identity 2x2 matrix. The symbol V» = 
denotes the spinorial covariant derivatives, being 
the spinorial connection Ej defined by Vil^ = 0. Then 
it leads to 



En = 0, 



e^ = \ht°t^ 



(8) 



Coming back to the Dirac equation ((3|), it takes the 
form 



r° f dt + °-h 



m = 0. 
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Restricting ourselves to the metric ((T|), the general so- 
lution of the latter equation consistent with Eq. ([2]) is 
given by 
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with arbitrary complex coefficients b±, 62, d\ and 
d,2- The only nonvanishing component of the energy- 
momentum tensor for the CDF is 



T D = ^(|6 1 | 2 + N 2 -N 2 -N 2 )-^, (11) 



where p DQ = a(b 2 - d 2 ), b 2 = |6i| 2 + \b 2 \ 2 and d 2 = 
\d\\ 2 + \d 2 \ 2 - For positive values of poo this source for- 
mally behaves as a perfect fluid representing a classical 
dust. However, the CDF will allow us to extend the 
analysis for negative values of the energy density. Here, 
we restrict ourselves to the physical sector d 2 < b 2 , this 
means that po > 0. 

From Eq. (|TT1) we see that the energy density of the 
CDF is given by p£> = ppo/a 3 . Thus we shall rewrite 
the Fricdmann equation ([2|) in the following form: 



3H 2 = p, n 



(12) 



Px 



PsO 



■ + A. 



(13) 
(14) 



Here a, b 2 and d 2 have dimensions of time -1 . 

Notice that there is some arbitrariness in the dis- 
tribution of the four components of Eq. pip between 
dark matter and dark energy, since it is not clear that 
we can represent each dark component with a mix of 
positive and negative energy components. To avoid 
this arbitrariness we have included both positive en- 
ergy components of CDF into dark matter and both 
negative energy components of CDF into dark energy 
This assignment has an invariant meaning because it is 
preserved under linear transformations of CDF. 

This distribution allows us to interpret the positive 
part of the Dirac energy-momentum tensor as a pres- 
sureless matter (giving rise to the total ("true") ob- 
servable matter p m ); while its negative part, together 
with the stiff fluid and the cosmological constant, we 
interpret as the effective dark energy component p x . It 
must be noted that this is equivalent to assume that 
we have a dust component pm = Paw a~ 3 and a CDF 
Pd = Pdo a -3 associated with a negative energy density 
(i.e. with d 2 > b 2 ) in the dark energy sector. 

The general solution of the Einstein equation (fT2|) 
with sources (Tl~3"|) and p4[) takes the form 



(10) a 3 (t) 




(15) 



where we have set the initial singularity at t = 0. In 
Fig [1] is shown the behavior of the scale factor and its 
derivatives. 



3 Dark energy evolution 

It is well known that the present observational data do 
not exclude the existence of the phantom dark energy; 
thus it is interesting to consider the possibility of cross- 
ing the phantom divide. We shall see that the Dirac 
field component can produce a strong variation of the 
dark energy density p Xl allowing one to have such a 
crossing at early times. 

In order to have a positive p x , we choose the parame- 
ters of the models according to the following restriction: 



2 d 4 < 4Ap s0 . 



(16) 
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Since a(t) is an increasing function of time, this effective 
dark energy component decreases until it reaches a min- 
imum value p xc = A — a 2 d 4 /4p s0 at a c = (2p s0 /ad 2 ) 1 ^ 3 
where the dark component crosses the phantom di- 
vide and begins to increase with time (see Fig. [2]). 
Fundamentally, the effective dark energy component 
crosses the phantom divide due to the presence of the 
d-parameter. In fact, the cosmological evolution of 
dark energy depends on the negative term —ad 2 /a 3 , see 
Eq. (|14p. because it produces a minimum at p x (a c ) = 0, 
showing the importance of considering the CDF as a 
source of the Einstein equation. 

Assuming that total matter and dark energy are cou- 
pled only gravitationally, then they are conserved sep- 
arately, so we have that 



pm + 3Hp m = 0, 
p x + 3H(l + w x )p x = 0, 



(17) 
(18) 



where we have assumed the equation of state p x — w x p x 
for the dark energy component. Taking into account 
that the dark energy component has a variable state 
parameter w x = w x (a), we define as a m = (p s o/A) 1//6 
the value of the scale factor where this equation of state 
coincides with the matter one, that is w x = 0. So that, 
the above restriction a 2 d A < 4Ap s0 now becomes a m < 
a c . In terms of the above parameters, a c and a m , the 
dark energy state parameter w x can be written as 



« ! 4 



1 - (a/a m ) e 



1 - 2(a/a c ) 3 + (a/o m ) e 



(19) 



On the other hand we can consider the state equation 
for the full source content. For this 2-component sys- 
tem we define the total pressure p T = w T p T , where 
P T = Px and the total energy density p T = p m + p x . 
This implies that 



UK 



1 



(20) 



where r = p m /p x . In Fig. [3] we compare the behaviors 
of the state parameters of the dark energy and the full 
source content. 

Let us consider in more detail the ratio of energy 
densities of matter and dark energy. The defined above 
ratio takes the form 



oh 1 a 



2„3 



Aa 6 - ad 2 c 



PsO 



(21) 



It is easy to show that this ratio has a maximum at 
a = a c = (pso/A) 1 / 6 (at this point the d 2 r/da 2 < 0). 
This maximum value is given by 



ab 2 y/p s0 /A 
2p s0 - ad 2 ^Psa/A 



(22) 




Fig. 1 We show the behavior of the scale factor (solid 
line) and its derivatives a (dotted line) and a (dashed line). 
From this it is clear that this cosmological scenario exhibits 
an accelerated expansion since there is a stage where a > 0. 




2 4 6 



Fig. 2 We show the behavior of p x as a function of 
time. The dashed line represents the limit CcLSC die — ££77-1, . 
i.e. a 2 d 4 = 4Ap s o, while solid and dotted lines repre- 
sent a typical case satisfying the condition a c > a m since 
a 2 d 4 < 4Ap s0 . 
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It is clear that for cosmological scenarios where r max < 
1 the dark energy component always dominates over the 
dark matter during all cosmological evolution. Thus, 
in order to have stages where dark matter dominates 
over dark energy, we have to require that r max > 1. 
So in this case we would have two values for the scale 
factor where the energy density of dark matter equals 
the energy density of the dark energy: 



a± 



a(b 2 + d 2 ) ± yja 2 {b 2 + d 2 ) 2 -AAp s0 



2A 



1/3 



(23) 



Thus, for cosmological scenarios where r max > 1, at 
the beginning the dark energy dominates over the dark 
matter until the scale factor reaches the value a = a_ 
where the dark matter energy density equals the energy 
density of dark energy and then it begins to dominate. 
This stage of domination of the dark matter is pro- 
longed until the moment when the scale factor reaches 
the value a = a + and the dark energy again starts to 
dominate over the dark matter (see Fig. |4|. 

3.1 Constraints on cosmological parameters 

In order to confront our model and the cosmological ob- 
servations we shall use the constraints on cosmological 
parameters obtained from the analysis of observational 
data assuming the A-CDM model. This choice is jus- 
tified by the fact that our model differs from A-CDM 
only at early times. It must be noted that in general 
such a procedure of using constraints derived from a fit 
of one specific model can give at best rough estimations 
for the parameters of a different cosmological model. 

The proposed scenario is characterized by four pa- 
rameters which may be constrained by the astrophys- 
ical observations available up to date. Since we have 
considered flat FRW cosmological scenarios, the dimcn- 
sionlcss density parameters are constrained today by 



"m,0 + "x,o — !• 

From Eqs. (TrJJ-fln]) we have that 

A. 



3ir = 



2 ab 2 : p s0 ad 2 



(24) 



(25) 



Evaluating it today (where we set a = 1) we have that 

p cHt = 3F 2 = ab 2 + p sQ - ad 2 + A, (26) 

so the two dimensionless density parameters are given 
by 



p m (a = 1) 

Pcrit 



ab 2 



ab 2 + p S Q — ad 2 + A ' 



(27) 



Px{a = 1) 



PsO — ad 2 + A 



(28) 



Pent ab 2 + pso - ad 2 + A ' 

Now it may be shown that, in general, for a flat FRW 
cosmology the deceleration parameter q is given by 



da 1 3 p 7 
~ = 2 + 2p~ 7 



a 



(29) 



Taking into account that we have p T = p x = u> x p x and 
p T = p m + p x the deceleration parameter is given by 



1 3 



(30) 



and evaluating it today (i.e. a = 1) and using (fl9)) we 
obtain 

a(b 2 - d 2 )(l - 2q ) + 2p sQ (2 - go) = 2A(1 + q ). (31) 

Thus, for accelerated scenarios, q < 0, we require a 
positive cosmological constant. Another constraint may 
be introduced by taking into account the moment when 
the Universe has started to accelerate again. In other 
words, this is related to the moment when the Universe 
starts violating the strong energy condition (p + p > 
and p + 3p > 0). So we must require the inequality 
p + 3p < 0. Now from the condition d = and the 
equivalent Friedmann equation 

£ = ~(p + 3p), (32) 
a 6 

we conclude that p + 3p = 0, which implies that p m + 
p x + 3uj x p x = 0, obtaining the condition 

a(b 2 - d 2 )(l + z acc f + 4p s0 (l + z acc f = 2A, (33) 

where the equation 1/a = (1 + z) was used. Here z acc 
is the value of the redshift when the Universe starts to 
accelerate again. 

In conclusion we have the four conditions ([2"T]) , (|2"8]l , 
(|3Tj) and ([3"3"| for the four parameters ab 2 , p s o, ad 2 and 
A of our model. 

Note that from Eqs. (fBTj) and (|33[) we have that 



(34) 



Pdo = a(b 2 - d 2 ) = Kp s0 , 

where 

4(1 + 9o )z 2 - 2(2 - go ) 



K 



(l-2 9 o)-(l + 9o)2 



, z = (i + z acc y c 



(35) 



Since p^o is related to the energy density of the CDF, 
we must require that K > 0. So from Eqs. (j27) . (j28]> . 
(I3T1) and (1331) we have that 



(36) 



ab 2 = 3-ffgO mj o, 



6 



PsO = 



w 2 



K(2 + z) + 2(1 + 2z 2 ) ' 



d 2 = b 2 - 



6KH$ 



t[K(2 + z) + 2(1 + 2z 2 )] 



A = 



3H 2 z{K + Az) 
A'(2 + z)+2(l + 2z 2 ) : 



(37) 



(38) 



(39) 



where we have used Eqs. (|24|) and (|34|). 

Now the four model parameters need to be con- 
strained. We do this by using the Eqs. and 
by considering the increasing bulk of observational data 
that have been accumulated during the past decade. 
The present expansion rate of the Universe is measured 
by the Hubble constant. From the fi nal results of the 



Hubb le Space Telescope Key Project (jFreedman et al 



2001) to measure the Hubble constant we know that its 
present value is constrained to be Hq = 68 ± 7 kms -1 
Mpc" 1 dChen et al. 112003ft . Equivalently we can write 



Hq 1 = 9.776 /i" 1 Gyr, where h is a dimensionless quan- 
tity and varies in the range 0.61 < h < 0.75. 

Now assuming a flat Universe, i.e. Eq. (I24p is valid, 



Perlmutter et al. ( Perlmutter et al. 1 119991 ) found that 



the dimensionless density parameter f2 mj o may be con- 
strained to be ~ 0.3, implying from Eg. ([24]) that 
Ur.n ~ 0-7 dCopeland et al. 1 120061: ICapozzielo 1 12007 ; 



Cardone et al. 



20041 ) . and the present day decelera- 
may be constrained to be —1 < 



<7o < —0.64 (Coueland et al. 


2006[ 


CaDOzzielo 


2007 


Cardone et al. 2004: Jackson 


2004 


• 



For consistency we also need to compare the age 
of the Universe determined from our model with 
the age of the oldest stellar populations, requiring 
that the Universe be older than these stellar popu- 
lations. Specifically, the age of th e Universe tp is 
constrained to be t p > 11-14 Gy r (ICopeland et al 



20061 iKrauss and Chabover I 120031 iGrattonet 2003 



1996tlRicher et al. Il200i lHansen et al. 112002ft . 

So let us calculate the age of the Universe from the 
Friedmann equation ([23]) . This equation may be writ- 
ten as 



H 2 



Hq ( Psa + aa 3 (b 2 



a 6 A 



p s0 + a(b 2 -d 2 ) + A 



(40) 



obtaining Eq. (|24|) when the expression (|40j) is evaluated 
today. Then the age of the Universe may be written as 



in 



r 

Jo 



dz 

77(T+I) 



1 r« 



a(b2- C P)+p s0 +A dz 
x(V-d*)(l+zf+p sa (l + z)»+A (1 + 2)' 



(41) 



In the Table [T] we include some values obtained from 
Eqs. (|36| (|39|) for the parameters ab 2 , ad 2 , p s o, A and 
K corresponding to some given values of the parame- 
ters Hq, qo and z acc (with c = 1 and G = 1). For the 
Hubble parameter Ho is considered both possible val- 
ues Hq- for h = 0.61 and Hq + for h = 0.75, and for the 
matter dimensionless density parameter we have taken 
^m,o = 0.3. The last two columns represent the age 
of the Universe determined from the model parameters 
for h = 0.75 and h = 0.61 respectively. Clearly in the 
proposed model there arc configurations which are al- 
lowed, being their ages larger than the oldest known 
stellar ages, since there exist combinations of the pa- 
rameters ab 2 , ad 2 , p s o, A which give to > 11-14 Gyr 
satisfying the stellar population constraints. 

3.2 The effect of the d-parameter 

Now it is interesting to get some insights concerning the 
nature of the d-parameter for studying the effect on the 
cosmological evolution of the negative term —ad 2 /a 3 in 



the Eq. (|14p for the dark energy component p x . It can 
be shown that in the proposed cosmological scenario the 
dominant energy condition (DEC) is violated thanks 
to the presence of this parameter. Effectively, if d = 
the dark energy state parameter w x and the state 
parameter of the full source content w T are given by 



PsO 



Aa 6 



PsO 



Aa 6 



PsO 



Aa 6 ' 



PsO 



Aa 6 + ab 2 a 3 ' 



(42) 



respectively. From these expressions we see that always 
— 1 < w x < 1 which implies that now the dark energy 
component satisfies DEC, as well as w T . Their general 
behavior is shown in Fig. [5] (compare with Fig. [3]). 

However, the fulfilment of the DEC does not imply 
that the Universe has a decelerated expansion. From 
Eq. ([32]) we may write that 



Imbriani et al. 112004 ; Hanssen et al. 12004 ; Jimenez et al. ^ 



65 
a 



-{P T + 3 Pt) 



Ap s0 , a(b 2 -d 2 ) 



-2A 



(43) 



and putting d = we see that the accelerated expansion 
is realized if a > a acc = ((ab 2 + V^ggg^o) )i/3. 
Another property of the d-parameter to be considered 
is its effect on the deceleration parameter qo. From 
Eq. (|38j) . which is independent of the Hubble parameter 
H and, using Eqs. (|3"5j) we can express the deceleration 
parameter go as a function of the z acc obtaining 



qo(z ac 



(z - 1) ((2z + l)(b 2 - d 2 )a - AH 2 (z + 1)) 
2iJ 2 (2z 2 + l) 



(44) 
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Table 1 In this table we show some values of the model parameters obtained for given Ho, qo and z acc (f2 m ,o = 0.3, c = 
1,G=1). 



Ho^ 1 ) 


qo 


Zacc 


cttf (s" 


■ 2 ) 


ctd A (s" 


■ 2 ) 




■ 2 ) 


A(s" 


2 ) 


K 


Ho+ , to 


Ho-, to 


Ho+ 


-0.68 


0.587 


5.315 x 10" 


36 


2.009 x 10" 


36 


2.368 x 10" 


37 


2.655 x 10" 


35 


14 


12 (Gyr) 


15 (Gyr) 


H - 


-0.68 


0.587 


3.516 x 10" 


-36 


1.329 x 10" 


-36 


1.566 x 10" 


37 


1.756 x 10" 


-35 


14 


12 (Gyr) 


15 (Gyr) 


H + 


-0.68 


0.94 


5.315 x 10" 


-36 


1.538 x 10" 


-36 


1.415 x 10" 


-39 


1.210 x 10" 


-34 


2619 


11 (Gyr) 


13 (Gyr) 


H - 


-0.68 


0.94 


3.516 x 10" 


36 


1.017 x 10" 


-36 


9.361 x 10" 


4(1 


8.006 x 10" 


-35 


2619 
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Fig. 3 We show the behavior of the dark energy state pa- 
rameter w x (dashed line) and the state parameter for the 
full source content w T (solid line). It is clear that the dark 
energy violates the dominant energy condition while the full 
source content does not violate it. Note also that the dark 
energy component remains in the phantom region as it en- 
ters into it after crossing the phantom divide. 




□ 



Fig. 4 We show the behavior of the ratio r = p m /px as 
a function of the scale factor a. The dashed line repre- 
sents the ACDM cosmology, where r oc 1/a 3 . The solid line 
represents the case where r > 1 and so at the beginning 
dark energy dominates, then at the range a £ (a_,a+) the 
dark matter component dominates over dark energy and, 
for a > a+, the dark energy component dominates again. 
The dotted line represents a scenario where the dark energy 
always dominates over the dark matter, since r ma x < 1- 
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where, as before, z = (1 + z acc ) 3 . We can see that in 
this case qo{z acc ) rapidly tends to the value 



So, 



2Hl 



<l> 2 



m.O 



b 2 



(45) 



obtaining, for d = and fi m> o = 0.3, the value qo,oo = 
—0.55 and, from this value the deceleration parameter 
reaches the value 90,00 = -1 for d w 6(see Fig.©. So 
the parameter d directly affects the range of validity 
of the deceleration parameter which is constrained to 
be in the range — 1 < qo < 0. Note that the value 
So, 00 = — 1 for d = b is independent of the value of the 
dimensionless density parameter fi mj o- 



4 Concluding remarks 

Since today the observations constrain the value of ui to 
be close to oj = — 1, we have considered broader cosmo- 
logical scenarios in which the equation of state of dark 
energy changes with time. The two principal ingredi- 
ents of the model arc a stiff fluid which dominates at 



early time ([Dutta and Scherrer 1120101 ) and a CDF. The 
positive part of the latter was associated with a dark 
matter component while its negative part was consid- 
ered as a part of the dark energy component and was 
responsible for the effective dark energy density cross- 
ing the phantom divide (see and compare Figs. |3] and [5J. 
At the end, this cosmological model becomes acceler- 
ated recovering the standard ACDM cosmology. 

In general, the model may be seen as a continua- 
tion of the inflation era. In the inflationary paradigm 
the scalar field driven by the potential V(4>), gen- 
erates the inflationary stage. In the slow roll limit 
$ 2 << V($), with w$ w —1, we have a superlumi- 
nal expansion while in the kinetic-energy dominated 
limit Q 2 >> V($>), with w$ ss 1, we have a stiff matter 
scenario characterized by a subluminal expansion. Tak- 
ing into account that our model has a variable equation 
of state, we can think of it as a transient model which 
interpolates smoothly between different barotropic eras 
as, for instance, radiation dominated era, matter domi- 
nated era an so on. In other words, from Eqs. (fT9|) -([2T |) 
we see that w x — > 1 (and w T — > 1) for a — > imply- 
ing that the energy density p x behaves like 1/a 6 and 
matching, after inflation, with the kinetic-energy mode 
of the scalar field p$ cx 1/a 6 . Now from Eq. (fT9|) wc 
see that p x passes through a radiation dominated stage 
(i.e. uj x = 1/3) for 



Q>rad 



2273 




(46) 




1-3 



Fig. 5 We show the behavior of w x and of w T as a func- 
tions of the scale factor. The dashed line represents the 
behavior of the dark energy state parameter, while the solid 
line represents the behavior of the state parameter of the 
full source content. Both satisfy the DEC since in this case 
d= 0. 



q 0.0 







1 2 

Zacc 






d = 
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d = 0.5 






d=1 





2.5 3.0 



Fig. 6 We show the behavior of the deceleration param- 
eter qo as a function of z acc for three values of the param- 
eter d (0, 1/2, 1) with Q m ,o = 0.3 and 6=1. We can see 
that the deceleration parameter rapidly tends to the values 
(-0.55, -0.6625, -1) respectively. 
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behaving like p x ss l/ a rad an d dominating over p m . 
After that, at a = a m , we have w x = and the effec- 
tive dark component behaves as a pressureless source 
thus obtaining a matter-dominated stage. Finally the 
model evolves from this state to a vacuum-energy dom- 
inated scenario. It is interesting to note that in the 
a matter-dominated stage, if the condition (fTBf is ful- 
filled, then the total energy density is given by p T = 
a(b 2 — d 2 )/a^ n + 2K and, if a, n = a c then p T = ab 2 /a^, 
being p x = 0, implying that we have at this stage only 
the dark matter component. 

The above results indicate that a cosmological sce- 
nario based on a CDF component and the effective 
multifluid configuration p x can, in certain cases, repro- 
duce the quintessential behavior (see Figs. [3] and [5]). In 
fact, the state parameter of the total matter content 
— 1 < w T < 1 is constrained the same as is the state 
parameter of the scalar field in quintessence models. In 
this manner we avoid the use of scalar fields and partic- 
ular classes of potentials for describing the dark energy 
component. 

Finally, all the parameters of the model have been 
expressed in terms of the observable quantities which 
may be constrained by the astrophysical observational 
data. In effect, in Table [1] some values of the model 
parameters ab 2 , ad 2 , p s o and A were included, which 
correspond to some given values of the parameters Hq, 
90i z acc and fl m ,o constrained by astrophysical observa- 
tions. 
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